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We propose an economical model of nonthermal leptogenesis following inflation during "instant" 
• preheating. The model involves only the inflaton field, the standard model Higgs, and the heavy 

' "right-handed" neutrino. 

o : 

^ , PACS numbers: 98.80.Cq 

o 

o 

■ I. INTRODUCTION 

' Leptogenesis is an attractive scenario to account for the observed matter-antimatter asymmetry of the universe. 
^ : In the scenario, a lepton asymmetry is generated by the decay of massive right-handed (Majorana) neutrinos, N, which 
are responsible for the (small) masses of left-handed neutrinos via the see-saw mechanism ^ . The lepton asymmetry 
Q^ ! is then translated to a baryon asymmetry by sphaleron processes around the electroweak era. Leptogenesis is 
' analogous to GUT baryogenesis 0| in that a massive particle decays with CP violation, and Boltzmann equations 
OO may be employed to track the time evolution of particles and the asymmetry 5J . The massive particles (in the case of 
leptogenesis, the A^s) must be created after inflation, either nonthermally or thermally during reheating, or thermally 
_ during the radiation-dominated era. Thermal leptogenesis 0, 0, 0, has been widely studied, but a careful analysis 
--^ , must be done in order to account for all the thermal and perturbative reactions (see Ref. Q and Ref. Nonthermal 

■ leptogenesis leptogenesis is an attractive alternative [lollll| . Here, we propose a new model of nonthermal leptogenesis 
involving instant preheating. 

Q ' Our model assumes hybrid inflation |l2j| . in which inflation is terminated by an abrupt transition in the properties 
^ ' of the scalar-field potential dominating the energy density during inflation. Such an abrupt change is often modeled 
^ as being triggered the action of a second "waterfall" field, causing the effective scalar-field to roll into another 
dimension in the scalar-field landscape. An important feature of hybrid inflation is that the properties of the scalar- 
field potential {e.g., the mass) during reheating may be quite different than the properties of the scalar-field potential 
during inflation.^ 

^ . , Our model assumes that the scalar-field energy is extracted and thermalized by instant preheating. In preheating 
^- El, particles are produced when the infiaton passes through a nonadiabatic phase around the minimum of the infiaton 
potential. In "instant" preheating Q , the inflaton is strongly coupled to a part icle whose mass depends on the value 
of the inflaton field. This particle can be either a boson 13] or a fcrmion 11 ; in our model we will assume it is a 
boson. As the inflaton oscillates, the coupling of the inflaton to the produced particle results in an increasing mass 
of the produced particle. As the mass of the produced particle increases, its decay rate will also increase, and decay 
channels disallowed when the produced particle is at the minimum of its potential may open. We will take advantage 
of both these properties. 

In our model we will assume that the inflaton couples to the standard-model Higgs boson, h, associated with 
electroweak symmetry breaking. We will also assume that, as expected, h couples to the N. Normally the mass of 
the Higgs, rrih, is much, much less than the mass of the N, ttin- However, during instant preheating this need not be 
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^ Although the scalar-field landscape may be quite complicated, involving several degrees of freedom, we will refer to the scalar field 
during preheating as the inflaton, although as remarked, the mass of the "inflaton" after inflation may be different than the mass of the 
inflaton during inflation. 



2 



the case, and the Higgs may decay directly into N, producing a lepton asymmetry. Later when the inflaton is close 
to its minimum, the produced A^s become heavier than the Higgs, and they will decay back to the Higgs. 

Let us elaborate on this picture. We will assume that the temperature, T, is always less than the mass of the N 
throughout the entire preheating stage. The mass of the Higgs will be determined by its coupling to the inflaton (/): 
nih oc \(f>\. For sufflciently large values of |0| during the inflaton oscillations, rrih will be larger than m^r. We will 
denote the absolute value of (f> when mh{(f>) — rriN 8S (j)c- 

It is useful to imagine a single oscillation of the inflaton field, in particular the first oscillation. As (/> passes near 
its minimum, h is effectively massless, and a burst of hs are created. The hs will decay to any kinematically allowed 
final states. Because of the large /i-top-quark coupling, the decay is predominately into top quarks, h ^ N becomes 
kinematically allowed when |0| becomes larger than 0^- Therefore, efficient lepton number production happens when 

is close to the minimum so h ^ N process takes place before all the h decays thermally into top quarks. In the case 
of hierarchical A^s, where nipf-^ = g\4>c\ ^ 'm-N2 — ffl02| ^ "mNs — ffl^sl, hs decay while |(/)| ^ \4>2,3\ due to the large 
/i-top-quark coupling. This process is nonthermal, as ruh > T at this time. Eventually (j) reaches a maximum point 
t/)™^'' and rolls back down. The decay of the h continues until (p < (pc- At this stage, N h decay happens, and hs 
continue to decay into fermions. A lepton asymmetry is generated by both h ^ N and N h decays. Another burst 
of hs are produced as </> passes again through the nonadiabatic phase at the origin, and the same events occur on the 
other side of the potential. Since h decays very rapidly, a negligible amount of hs remain when (f> repasses through 
the nonadiabatic regime to produce more hs. This eliminates the influence of the old hs with during production of 
new hs, and the backreaction of Higgs in the nonadiabatic region need not be considered. The production and decay 
of hs siphon away energy from ip, and (p'^^^ decreases for each oscillation. A schematic diagram of the regions of the 
potential in instant preheating is shown in Fig. ^ 

There have been many models of leptogenesis. A hallmark of our model is the economy of fields. The only 
undiscovered fields are the inflaton, cp, the standard model Higgs, h, and the right-handed neutrino, N. There are 
very good reasons for suspecting that all exist! The only unfamiliar aspect of our model is the strong coupling of the 
inflaton field to the Higgs field. While there is no reason to preclude such a coupling, it would be very interesting to 
find particle-physics models with a motivation for the coupling. 

In the next section we present the Lagrangian used in our calculation, we parameterize the nonadiabatic creation 
of particles in preheating, and we discuss the decay rates and CP violation parameters. In Section HI we present and 
solve the Boltzmann equations used in the calculation, presenting the main results in the form of figures. The final 
section contains our conclusions. 



II. INSTANT PREHEATING AND THE SEE-SAW MECHANISM 



Inflation ends in the hybrid model when </> meets a "waterfall" potential in another direction of the scalar field 
landscape. The (p promptly falls into this potential which is responsible for preheating. Hence, (p does not carry 
restrictions on potential parameters (such as the mass) deduced from present cosmological observations. For instance, 
the mass of the (p during the preheating process may be more massive than the mass of the (p during inflation. 

We assume cp is coupled to the standard model Higgs h, with interaction Lagrangian of the preheat field given 
by 

^S^'prcheat = -^ff'^'/^', (1) 

where g is the coupling constant. Ignoring its electroweak-scale mass, ruh = g\(p\- We define 4>c = g/fnN- Thus, 
depending on the initial condition of the (p field, nih may become larger or smaller than tun as <p oscillates about the 
minimum of its potential. 

The inflaton-Higgs coupling leads to a potential of (p about the minimum in the form 

m = i/iV + ^ffV/^^ (2) 

where is the (p mass. The Hartree approximation will be used to take the average value of jl^ . where we 
approximate « {h?)- Of course h^ is formally inflnite, but becomes 2nh/mh after renormalization. The equation 
of motion of (p becomes 

4> + 3H4> + fi^(p + 2gnh(p/\(p\=0, (3) 
where the dot stands for the time derivative, H = d/a is the Hubble expansion parameter, and a is the scale factor. 
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FIG. 1: A schematic diagram of regions in the infiaton potential during instant preheating. The shaded column around the 
minimum illustrates the nonadiabatic region where hs are created. In regions of \(l>\ > (l>c, h ^ N decay occurs. In regions of 
101 < <l)c (modulo the nonadiabatic region), N h decay occurs. The regions are not drawn to scale. 



The hs are created when </> goes through a nonadiabatic phase, which occurs near the minimum of the potential [13| . 
This phase is very short and can be treated as instantaneous. The number density of hs created in the nonadiabatic 
phase is HI [13 

„,.,0) = ^ f .M"„. ^^^f-lk = (4) 

where (jjQ is the initial time derivative of (f), and k is the momentum. A large coupling constant g 1 enables a quick 
and effective thermalization of the universe within a few oscillations of 4>. 

The see-saw mechanism Lagrangian for three families with Majorana neutrino masses mpf. (i = 1, 2, 3) and Yukawa 
couplings to the Higgs and light neutrinos I is given by 

i?see-saw = Nf + Y^j^^UN^h. (5) 

The left-handed light neutrino masses are nii, — —{vY^)'^ni^{vY^)^ where v — 247 GeV is the Higgs vacuum 
expectation value, ^sec-saw also generates a dimension-5 effective operator which causes CP violation among the 
leptons. Throughout the paper we consider the case of hierarchical Majorana neutrinos, m^Vi <SC 'm^^ <^ niN^. Hence 
the CP asymmetry is created during {(pd < \<P\ \4'2,3\, as the h Nl interaction competes against the dominant 
background interaction of /i — s- //. This mass hierarchy allows us to consider only interactions involving Ni; hence 
we will drop the family subscript unless distinction between A'^i and 7V2,3 is required. 
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The decay processes h — > Nl{l) or N ^ hl{l) give rise to a lepton asymmetry. Both processes are possible in our 
model: the former process happening when nih > mjv {i.e., when \(j)\ > (j}^), and the latter happening for mh < tojv 
( i.e., for \(f>\ < (j)c). The CP parameters in these interactions, Ch and cn, respectively, are defined as 



eh 



h^Nl 



h^Nl , 



r 



r 



N^hl 



N^hl 



h~*Nl 



h^Nl 



^N^hi + r 



(6) 



N^hl 



where the subscripts of the decay width F denotes the decay process concerned. The possible combination of these 
CP violating processes are 



Nl ■ 



Nl 



Ml 
hll 
hll 
hll 



(7) 



The second and third final states have zero lepton number, but the first and fourth final states give rise to lepton 
asymmetry. The total CP asymmetry etot is expressed as 



Etot 



h^Nl 



r 



h^Nl 



N^hl 



h^Nl 



h^Nl 



I ^N^hi + r 



The explicit expression of this CP parameter e at one-loop is (1| 



-2E 



(rtr)ii 2P-Pi 



(8) 



(9) 



for both h and A'' decay, /q and /i are the tree- and one-loop diagrams, respectively. In Eq. ®, P is the four- 
momentum for either h or N , and Pi is the four-momentum of I. One of the condition of our leptogenesis model 
is the nonthermal production of N . All processes involving CP violation are nonthermal and a zero-temperature 
expression of e is required. Because of the mass hierarchy, A^2,3 are not on-shell in the one-loop diagram. The h and 
I are considered massless {mh^i ^ tijv) in calculating ejv, and N and / are considered massless (m^v,; <C rrih) for the 
calculation of eh. Therefore the expression of the above equation is equivalent for both e/i and e^r, and thus — tN. 
We look into the case oi tN for the N ^ h decay in the following as an example. 

In a hierarchical iV- family structure (m^Vi ^ ™JV2 3)1 the zero-temperature expression of Eq. Q becomes 



im[(ytr)2^ 



''Ni 



/(^) 



A convenient parameter to use is the effective neutrino mass 



(1 +a;)ln 



1 



(10) 



ruN 



(11) 



which is the contribution to the neutrino mass mediated by A'^i. A way to understand mi is to use an example of a 
hierarchical lefthanded neutrino spectrum of mi <C m2 = mgun ^ "Za = matm where msun and matm are the deduced 
solar and atmospheric neutrino mass from neutrino oscillations. If A^i gives rise to the atmospheric mass splitting 
then mi = matm; if Ni causes the solar mass splitting then mi > msun- 



The explicit expression of ejv is [Tsl l 



1^3 rriN [ms - mi) 
< 5 X 

16 TT 



1 — mi/fhi if mi <^ 
^1 — m\/rh\ if mi ~ 



ma, has not been measured. We assume ma — max(mi, matm) = 0.05 eV. Eq. H12(l is maximal when mi = 0; 

I max, 3 ITLN^ m3 



16 TT 



(12) 



(13) 



The above expression holds for |e™'''^| as well. 
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The h decay thermalizes the universe. Before electroweak symmetry breaking in the standard model, only the 
fermion decay channel is allowed. The decay width of ft. — *■ // is given by 

9 

r.^// = 2^^"'^' (14) 

where m/ is the mass of the fermion. The decay widths of h Nl and h — s- Nl differs only by the CP asymmetry. 
When calculating the Boltzmann equations, the CP asymmetry is factored out as e and an identical expression for 
^h^Ni E^nd r^^^;- are used. In these cases, II will be dropped from the reaction expressions. The same holds for 
N ^ hi and hi. The decay widths for these are 

r;,_^ = iY^Y)n'^, (15) 
F^^. = {Y^Y)n^. (16) 



III. LEPTOGENESIS 



We study the time evolution of the hs, Ns, and the lepton asymmetry by means of the Boltzmann equations. Along 
with the equation of motion for in Eq. |j2Jl, the following set of Boltzmann equations are used: 

nh + iHnh + Tf^_^fj:{ni^ - nl^) + rh^N{nh - nl'') -TN^hiriN ~ n''^) = 0, (17) 
UN + iHnN + TN~.h{nN ~n']^) ~Th^Ninh-nl^) = 0, (18) 

ul + SHriL - Y^h^N{nh ~ n^) - -^^N^h{nN - "-^) = , (19) 

PR + AHpR~Tf^^j:f{nh-nl^)-Th^Nm,i{nh~nl^)-TN^h{nN-n''j^)rnN = 0. (20) 

Here, njv is the number density of N, = ni — nj is the lepton number density, and pi^ is the radiation energy 
density. It is understood that T^^n occurs when to^ > m^r, and Fjv— ►/! occurs when m/i < niN. The expansion rate 
H is 



H' = ^^(^l4>^ + V{c^)+Ph + PN + PR], (21) 



where Mpi is the Planck mass, and ph and pN the h and N energy densities density. 

There is an implicit assumption of a rapid thermalization from and h decay products. The h decay products 
include fermions (Eq. 114|) '). where the dominant channel are decay into top quarks, which have large cross section and 
will result in rapid thermalization. Indeed, instant preheating relies on rapid thermalization from the large coupling 
of <j)^-h^ and fast decay of hs [T^ . 

The AL = 2 off-shell scattering is not included in Eq. (|19|) . but as its decay rate is very small compared to the 
on-shell scattering Q, we neglect this effect. 

The hs are partially thermal throughout preheating. They are nonthermal when |0| > 0c, i-e., rrih > rriN- They 
may be thermal when < <f>c, depending on the thermalization rate, and its effect must be taken into account. The 
equilibrium number density of h, using Maxwell-Boltzmann statistics, is 



27r V T 



K^imh/T), (22) 



where K2 is the modified Bessel function of the second kind. The Ns are nonthermal throughout preheating. Conse- 
quently S> ?T.jv, and terms are neglected in the above Boltzmann equations. 

We now express Eqs. Q and H17|I - I|2U|) in terms of dimensionless variables and convert time derivatives to derivatives 
with respect to a; = pt. The dimensionless variables are 

N^^ ^ ^A'; Nu^^A^, Nr^^^A'- N,^^A\ (23) 
nh{Q) n/i(0) 7i/i(0) n,i(0) 
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Here, hq is the initial scale factor. The dimensionless equation of motion of cj) and the Boltzmann equations become 

^" + 3W + ^ + p(^^~^^ = 0, (24) 

K - {^K - K') + {Nh - K") - ^^^Nm = 0, (25) 

^ /I 

„ (iV, - iV,f ) + E^Nm = 0, (26) 

N'L-'i {Nh - K') ~ '-^^^Nn = 0, (27) 

R' _ £^^^(^, _ N'^'>)A - El^(3^{Nn - K'')A - E^^I^NnA = 0, (28) 
H fj, fj, 47r (jj^' 

A' -An ^ 0, (29) 

where the prime superscript denotes derivative with respect to x, and /3 = g^^'^ (j)]/"^ / {Ai:^ . 

The range of value used for the parameters during the numerical integration are as follows: 3 x 10~^ eV < rhi < 1 
eV; 10^ GeV < tojv < lO^^ GeV; fi > lO^^ GeV; < (10^^ GeV)^; and g ~ 1. The upper limit to comes from the 
sum of the three left-handed neutrino mass combining neutrino oscillation data with constraints from the cosmic 
microwave background and large scale structure observations 16,], under the assumption that rhi < ^rrii, with a 
hierarchical left-handed neutrino spectrum. The lower limit has been arbitrarily set. The Yukawa coupling must be 
neither too small nor too large for the see-saw mechanism to be compelling. The upper bound of mAr is derived 
from Eq. by setting {YjY,y)ii ~ 100. The inflation parameters /i and 4'o are derived from observation In 
the preheating model we consider, the mass of the inflaton during preheating must be larger or equal to the mass of 
the inflaton during inflation. Hence we consider fi > 10^^ GeV. We stress that the bounds of all the parameters are 
approximate and not very stringent. 

We terminate preheating when pr/ > 10, deeming this to be sufhcient that the radiation energy dominates over 
the scalar energy density. An example of numerical results for a model is shown in Fig. |21 

Reheating is completed after twelve oscillations for this case, but only the first few oscillations are shown for the 
sake of clarity. It is clearly seen that h decays very quickly after it is created. Most of the decay happens in the 
h ^ ff channel. For efficient leptogenesis, (pc must be close to the origin in order for the h ^ N process to start 
before the hs all decay away. Because of the quick decay of the /is, each oscillation can be treated as independent from 
each other. The decrease of the maximum point of Lp indicates the energy of the (j) field going into thermalization of 
the universe. is normalized with e{Y'^Y) in the figure. 

We have calculated the lepton number n^/s , where s is the entropy. The lepton number gets translated into a 
baryon number ub/s by sphalcron processes. Sphalerons transfer a lepton asymmetry to a baryon asymmetry by 
reactions conserving ub-l but violating ub+l- The relation between baryon number and lepton number is ,19] 

ns _ ^ riB-L _ C riL 24 + 45,, 

s s G — Is 66 4- log/i 

where gu is the number of h generations. We consider a one -gen eration h model. The observation of baryon number 
comes from big bang nucleosynthesis (BBN) considerations |20j |. cosmic microwave background determinations |2lj |. 
and large scale structure measurements observations, with an assumed cosmological model. In a ACDM cosmology, 
these observations imply 

!^ = 6.5+Oi X 10-1° ^ ^ «9xl0-" (s = 7.04n^). (31) 

Figure |21 shows the region of niN and mi where ns / s is higher than observation. We have used the maximally 
allowed CP parameter |e""*^| in our Boltzmann equations, but |e| does not alway retain the maximal value. The lower 
limit of rhi is due to the bound of 3 x 10^^ eV we used in our calculations; if the bound is lowered, the contour simply 
continues downward. The slant shape on the left hand boarder of the shaded area is not a simple slope relation; this 
comes from the combined restriction of Hb/s > 9 x 10"^^ and the nonthermal condition of mjv > T. The preheat 
field parameters /i and 0o are not very sensitive in determining ub/s. 

The reheat temperature Trh is greater than lO^" GeV in most of these regions. In supersymmetric models, this 
leads to overproduction of gravitinos which causes incompatibility with BBN observations 23]. Some models of 
supersymmetry have a larger mass to the gravitino |24| . which can relax the constraint on Trb- As we do not 
explicitly consider supersymmetry in our calculations above, our model agrees with all observations. 
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FIG. 2: Evolution of Kl = NL/{e{Y''Y)ii) (solid line) and R (dashed line) in the upper box, \(p\ (dashed curve) and Nh (solid 
line) in the lower box, during preheating. Most of the lepton number is created during the early stage of oscillations. Note 
that h decays very quickly. Parameters used are 00 = 10 GeV , fj, — 10 GeV, mjv — 10 GeV, mi = 10-" eV, and 5 = 1. 



IV. CONCLUSION 



We have proposed a simple, economical model of nonthermal leptogenesis during instant preheating in the context 
of standard model and its extension to include Majorana partners. A hybrid inflation is employed, which allows us 
to evade the constraints on the properties of the inflaton potential from observations. 

The assumed strong coupling between and h ensures a quick thermalization. Preheating occurs within a few 
oscillations justifying the definition of "instant" in this preheating scenario. The dominant thermalization process is 
h — > //. In addition to this process, the see-saw mechanism provides another decay channel via the N-h coupling. As 
nifi cx 0, rrih can be larger or smaller than niN throughout preheating. When rn/j > ttin^ h N decay occurs, and the 
opposite happens when irih < tojv- Both processes produce lepton number. The lepton asymmetry is produced under 
nonthermal conditions, since the T < tojv- The lepton number subsequently gets transformed to baryon number via 
sphaleron process. The effective light neutrino mass rhi has an upper bound of 1 eV. For a successful leptogenesis to 
happen in our model, we require mpf > 10^^ Gev. For most of the parameter space we find Trh > 10^'' GeV, which 
may cause incompatibility with BBN observations in some SUSY models. 

To summarize, if the electroweak Higgs is coupled to the inflaton, then one can expect instant preheating where 
the inflaton energy is extracted by resonant Higgs production as the inflaton passes through = 0. As the inflaton 
grows during an oscillation, the effective mass of the Higgs may become large enough such that it can decay to the 
right-handed Majorana neutrino N, even if the mass of the N is as large as lO^'^ to 10^^ GeV. A lepton number may 
be produced in this phase. Later, when the value of the inflaton fleld decreases, the Higgs mass will fall below the N 
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FIG. 3: Regions of mjv and rhi which satisfies the observed ns/s value. 

mass, and the N will decay to Higgs, also producing a lepton number. 

The resulting lepton number only weakly depends on inflation parameters, is rather more sensitive to two neutrino 
mass parameters from the neutrino sector, and depends on the CP-phase in the heavy neutrino sector. 

Our model is yet another scenario for leptogenesis, and illustrates the cosmological richness of the well motivated 
see-saw explanation for neutrino oscillations. 
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